Axially and spherically symmetric solitons in warm plasma 
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We study the existence of stable axially and spherically symmetric plasma structures on the basis 
of the new nonlinear Schrodinger equation (NLSE) accounting for nonlocal electron nonlinearities. 
The numerical solutions of NLSE having the form of spatial solitions are obtained and their stability 
is analyzed. We discuss the possible application of the obtained results to the theoretical description 
of natural plasmoids in the atmosphere. 
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I. INTRODUCTION 

The studies of stable spatial solitons is an important 
problem of contemporary physics There are numer- 
ous manifestations of stable solitonic solutions of non- 
linear equations in two and three spatial dimensions in 
nonlinear optics [1] , solid state Q and plasma physics [1] . 

In plasma physics it was established [5j that the nonlin- 
ear electron-ion interaction leads to the modulation insta- 
bility and results in the collapse of a Langmiur wave . 
In contrast to the electron-ion interaction, the nonlin- 
ear electron-electron interactions, studied in Refs. @, @], 
were shown to stabilize the evolution of a Langmuir wave 
packet, making possible the existence of stable spatial 
plasma structures. 

It is convenient to describe the evolution of nonlin- 
ear Langmuir wave packets on the basis of a nonlinear 
Schrodinger equation (NLSE) The nonlinear terms 
studied in Refs. @, Q are local since their contributions 
to NLSE contain only a certain power of the electric 
field amplitude. The nonlocal terms in NLSE, derived 
in Refs. (ToL [Tl| . are also important when, e.g., a wave 
packet is very steep. Under certain conditions these non- 
local electron nonlinearities can arrest the Lamgmuir col- 
lapse. It is suggested in Ref. [ill that a nonlocal NLSE is 
a theoretical model for stable spatial plasma structures 
obtained in a laboratory [l^ . 

Besides the nonlocal electron nonlinearities taken into 
account in Ref. [ll| there are analogous contributions to 
NLSE originating from the electron pressure term. These 
terms have the same order of magnitude as those con- 
sidered in Ref. [ll] and are important for plasma with 
nonzero electron temperature. Note that nonlinear waves 
in warm plasma were also studied in Ref. |13| . In the 
present work we carefully study these additional nonlin- 
earities and examine their contribution to NLSE. 

This paper is organized as follows. In Sec. [Til on the 
basis of the system of nonlinear plasma equations, we ex- 



amine electrostatic plasma waves having axial and spher- 
ical symmetry. Then, in Sec. Ill Al we briefly review the 
previous studies of the influence of electron nonlinearities 
on the dynamics of Langmuir waves in plasma. The new 
NLSE, describing stable spatial plasma structures, tak- 
ing into account the nonlocal nonlinear terms is derived 
in Sec. Ill Bl In Sec. Mil we analyze solutions of this NLSE 
numerically. We consider the possible application of our 
results to the theoretical explanation of the existence of 
atmospheric and ionospheric plasmoids in Sec. IIVI Fi- 
nally we summarize our results in Sec. [V] 

The description of plasma waves in frames of Lagrange 
variables is presented in Appendix [X] 



II. THE DYNAMICS OF SPATIAL LANGMUIR 
SOLITONS 

In this section we study axially and spherically sym- 
metric waves in plasma accounting for local and nonlocal 
electron nonlinearities. We derive a new NLSE which is 
shown to have solitonic solutions. 

To describe electrostatic waves in isotropic warm 
plasma, in which the magnetic field is equal to zero, 
B = 0, we start from the system of nonlinear hydro- 
dynamic equations, 
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V • (n e v e ) = 0, 
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where the densities of electrons and ions, v e ,j are 

their velocities, E is the amplitude of the electric field, 
m is the electron mass, e > is the proton charge, and 
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is the pressure tensor which is calculated using the elec- 
tron distribution function f e . Note that Eq. (|2.1|) follows 
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from more general Vlasov kinetic equation for the func- 
tion f e . 

It is known that Eq. (|2.ip allows small amplitude Lang- 
muir oscillations on the frequency lj p = y 4irn e 2 /m, 
where no is the unperturbed electron density. However, 
since Eq. (|2.ip is nonlinear, the higher harmonics gener- 
ation is possible. Therefore we can look for the solution 
of Eq. (|2.ip in the following form: 
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amplitude Ei was obtained: 
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where we separate different time scales, 
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The functions n s , E s , and v s in Eq. (2.3) and the ampli- 
tude functions 711,2, Vi,2 and Ei, 2 in Eq. (|2.4|) are sup- 
posed to vary slowly on the l/(J p time scale. Moreover 
we suggest that, e.g., no 3> «i 3> n s ,2 etc, i.e. slowly 
varying functions, marked with index "s" , and the am- 
plitudes of the second harmonic are much smaller than 
the corresponding amplitudes of the main oscillation. 

We neglect the velocity of ions in Eq. (|2.ip and suggest 
that ion density is represented as n^ = no + n, where the 
perturbation n is also a slowly varying function on the 
1/ojp time scale. Note that n does not necessarily coincide 
with n s . 

We will study electrostatic plasma oscillations in two 
or three dimensions. Thus we assume radially symmetric 
quantities in Eq. (|2.ip . 



= e r x 



v e (r,t) 
E(r,t) 



_ ( n e (r,t) 
m(r,t) 



, (2.5) 



where r is the radial coordinate and e r is the basis vector 
in spherical or cylindrical coordinate system. 



A. Local electron nonlinearities 

The contribution of electron nonlinearities to the evo- 
lution of Langmiur waves in frames of the model (|2.ip ~ 
(|2.5p was taken into account in Refs. [7], [f| and the fol- 
lowing equation for the description of the main oscillation 



iEr+^rfjVCV • Ex) - 
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(3(d) Ei|Ei| 
127rmnoW p r 2 



= 0, 



(2.6) 



where td = ^/T e /47re 2 no is the Debye length, T e is the 
electron temperature, and j3(d) = (d — 1)(4 — d)/2, d is 
the dimension of space. 

Eq. (|2.6p should be supplied with the wave equation 
for the ion motion 14], 
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(2.7) 



where c s = y (T e + jjT)J/M is the sound velocity, Tj is 
the ions temperature, 7, is the heat capacity ratio for 
ions, M is the ion mass, and A is the Laplace operator. 

At the absence of electron nonlinearities [the last term 
in Eq. (|2.6p ] the system (12. 6p and (|2.7[) corresponds to 
the Zakharov equations [5( . It should be noted that the 
contribution of local electron nonlinearities is washed out 
from Eq. (|2.6I) in one dimensional case d = 1. 

The Zakharov equations are known to reveal the col- 
lapse of a Langmuir wave packet @ : the size of a wave 
packet is contracting and the amplitude of the electric 
field is growing. It was shown in Refs. @, @ that Lang- 
muir collapse can be arrested and stable spatial plasma 
structures can appear in two and three dimensions, d — 
2, 3, since the second nonlinear term in Eq. (|2.6p is defo- 
cusing. 



B. Nonlocal electron nonlinearities 

The influence of electron nonlinearities on the dynam- 
ics of a Langmuir collapse was further studied in Ref . [ll[ . 
Using the relation between slowly varying electron den- 
sity n s and the perturbation of ion density n, 



AjEij 

4ne 2 



+ f D An s 



(2.8) 



which was obtained in Ref. [8j , we can approximately find 

n s as 



AjEij 

47TC 2 



(2.9) 



Note that the last term in Eq. (|2.9p . ~ r^An, which is 
important at rapidly varying ion density, was omitted in 
Refs. 

Using Eqs. ([2~7]) . (|2~9| and supposing that T t < T e , 
one obtains the new nonlinear term in the left hand side 
of Eq. (EU) (see Ref. [HI), 



EiA|E!| 2 



(2.10) 
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This new contribution was shown in Ref. [11| to arrest 
the Langmuir collapse. Moreover the nonlocal nonlinear- 
ity (|2.10j) is more effective in preventing the collapse com- 
pared to that found in Rcfs. [7, 8]. It should be also noted 
that the new nonlinear term predicted in Ref. does 
not disappear in one dimensional case. The nonlinearity 
analogous to that in Eq. (I2.10j) . ~ EiA|Ei| 2 , appears in 
the Zakharov equations with quantum effects [15[ , while 
taking into account the quantum Bohm potential. 

The nonlocal electron nonlinearities, analogous to that 
studied in Ref. [11], can follow not only from Eq. (|2.9|) . 
We can consider the contribution of the slowly varying 
electron density to the electron pressure (|2.2I) . Analo- 
gously to Eqs. (|2.3[) and (|2.4p one can discuss the de- 
composition of the distribution function f e proposed in 
Ref. 0, 



fe =/o + fs + he 



2itu v t , 



(2.11) 



where fo is the equilibrium distribution function. For 
classical plasma it can be, e.g., a Maxwell distribution 
corresponding to T e . All the quantities in the expansion 
series (|2.11[) depend on v — v 



(i) 
/ 



Using the result of Ref. |7| we represent /i as, 

i \df df s 
Ji= — -K- v ' Vjvi + — (v • Vjvi 

-(vi-V)/ s + -E s ^i + 
m av 



(2.12) 



where we drop terms containing fa, E2, and higher power 
of the first harmonic amplitudes. The function f s in 
Eq. (|2.12[) was also found in Ref. 7] for the case of 
isotropic equilibrium distribution, 



-v T n 6 



Idfo 
v dv 



(2.13) 



where i>t = y/T e /m is the thermal velocity of electrons. 
The terms which do not contain the slowly varying den- 
sity n s are omitted in Eq. (|2.13[) . 

Now we can express the nonlinear term \7iPij/n e in 
Eq. (EU as 
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Using the Maxwell equilibrium distribution function, 



fo = «o 



\2nT e 



3/2 



exp 



rav 
27; 



(2.15) 



normalized on the unperturbed electron density, we can 
express the gradient of pressure in Eq. (|2.14l) in the fol- 



lowing form: 



cv. 



^-e~ i ^ t [W 3 (n s W l E u ) +V i (n s ViE lj ) 

+ V l (n s V 3 E li ) + Vj(E u Vin s ) 

- {Vjn s )(ViE u ) - 3n s Vj(ViE u )], (2.16) 

where we use the relations between the amplitudes of the 
main harmonic, 



Vi 



-Ei (V • Ei) = -47rem, 



(2.17) 



established in Ref. M- The leading term in Eq. (|2.16[) 
was derived in Refs. [3, H|. The next-to- leading terms, 
proportional to the derivatives of n s , are important when 
one has rapidly varying in space wave packets. 

Using Eqs. (|2.9[) and (|2.16|) we can obtain the gener- 
alization of Eq. (|2.6p which takes into account the non- 
local electron nonlineariries due to the interaction of a 
Langmuir wave with the low-frequency perturbation of 
electron density, 

• 3 , , w„ B(d) EIEI 2 

2 2n lzirmnoUJp r z 
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■ [EAn - V[n(V • E)] - V 2 (nV 4 E) 



2noOJ p 

- Vi(riVEi) - V[(E ■ V)n] + (V • E)Vn 

+ 3nV(V ■ E)] = 0, (2.18) 

where for simplicity we omit the index "1": E = Ei. 
Eq. (|2.18j) should be supplied with Eq. (|2.7j) governing 
the evolution of ion density perturbation n. 

The nonlocal term ~ EAn in Eq. (|2.18l) was derived in 



Ref. 11(. The remaining nonlocal nonlinearities, which 
are of the same order of magnitude as the term ~ EAn, 
were omitted in that work. 

Let us study the evolution of the system (|2.7[) 
and (|2.18[) in the subsonic regime, when one can neglect 
the second time derivative of the ion density in Eq. (12.71) . 
Considering axially symmetric, d = 2, or spherically sym- 
metric, d — 3, cases [see Eq. (|2.5|) ] and using the dimen- 
sionless variables, 



3 r 9 E 

- — , s = -(j p t, $ = =, (2.19) 

2r D ' 4 p ^18Trn Q T e V ' 



we can represent the dynamics of the system (|2.7[) 
and (|2.18[) in the form of a single NLSE, 
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It should be noted that Eqs. (j2~T8)) or (j2~20)) does not 
conserve the number of plasmons, 



(2.21) 



A = / dV\E\*, A ^ 0. 



d^/dx 



This fact is because of the presence of the term 

in Eq. flUO]) . 

To analyze the integrals of Eq. (|2.20[) let us separate 
the variables in Eq. (|2.20|) : $ = e~ iXs d>(x). Then mak- 
ing the following nonlinear gauge transformation of the 
"nonhermitian" Eq. (|2~2D|) : 



xjj = cxp 



-> 2 



we can cast it in the equivalent form, 
1 d 



d_ 

dx 
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(2.22) 
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(2.23) 



where = d 2 /dx 2 + (d-l)/ix d/dx is the radial part 
of the Laplace operator. To derive Eq. (|2.23[) we use 
the decomposition of the transformed "wave function" 
squared Q, |Vf - (1 - ||0| 2 + • • ■ )|0| 2 » |</f, and 
keep only cubic terms in Eq. (|2.23l) . 

The modified NLSE (f2T23|) is "hermitian" , i.e. it con- 
serves the number of plasmons, 



N,i, = 0„ 



da; x 1 



d-l 



(2.24) 



where we restore the s dependence, ^ 



-i\s 



ip(x), and 

ild = 27r, for d = 2, or Qd — 47r, for d — 3, is the solid 
angle. It should be also noticed that the nonlocal non- 
linearities does not disappear in one dimensional case as 
the local ones studied in Rcf. 0, Q. The dynamics of 
Langmuir waves in warm plasma in arbitrary dimensions 
is analyzed in Appendix using Lagrange variables. It 
is shown there that in case of nonzero electron temper- 
ature the nonlinear terms does not disappear in higher 
dimensions d > 1. 

We can notice that Eq. (12.23[) is analogous to NLSE 
equation derived in Ref. It also contains ipA x \i/;\ 2 

term, although the coefficient is different. The main 
discrepancy is the presence of the term ~ (d — l)/x x 
ipd\ip\ 2 /dx. In Sec. IIIII we analyze the influence of this 
new contribution numerically. 



III. NUMERICAL SIMULATION 

It is difficult to construct other conserved integrals of 
Eq. (|2.23[) . e.g., a Hamiltonian, independent of the num- 
ber of plasmons (|2.24|) . Therefore one has to analyze the 
behaviour of solutions of this equation numerically. 



Eq. (|2.23[) should be supplied with the boundary condi- 
tions, ip(0) — ip(oo) — 0, and thus treated as a boundary 
condition problem. We have found numerical solutions of 
this problem using a boundary condition problem solver, 
incorporated in the MATLAB 7.6 program. It requires 
an initial "guess" function which was chosen as 



ip g (x) = y/NoA d xexp 



X 

2^2 



(3.1) 



where a is the "width" of the function, Ad — 1/(\/tt<t 2 ), 
for d = 2, and A d = y/2/3/ \n 3 / 4 a^ 2 ), for d = 3. The 
function (|3.1[) is now normalized on the initial number of 
plasmons Nq, to be compares with the actual number of 
plasmons Nj, obtained from a numerical solution. 

The best convergence of the numerical procedure is 
achieved when a corresponds to the minimal value of A 
in Eq. (|2.23l) at the given Nq. This analysis is analogous 
to the trial function method [l6j for minimizing of the 
Hamiltonian of NLSE. In Fig. |TJ(a), for d = 2, and in 
Fig. [2Ja), for d = 3, we present the dependence of A 
versus a for different values of Ao- One can see that 
the function A(<r) has a minimum if N > A cr . In two 
dimensional case the critical number of plasmons can be 
easily found: N CI = 8tt. Note that one can expect the 
stability of a soliton with respect to the collapse if A < 
for d = 2, whereas in 3D case there is some range of 
positive A which corresponds to uncollapsing solutions. 

The solutions which correspond to various values of A, 
a, and are shown in Fig. QJc,d), for d = 2, and in 
Fig. [2Ic,d), for d = 3. A "guess" function which does 
not correspond to the minimum of the function A(cr) also 
gives some solitonic solution of Eq. (|2.23p . however the 
convergence is much worse than in the minimal A case. 
Indeed, significantly differs from Ao for such a so- 
lution. When the deviation from the minimal A is big, 
although A remains to be negative, no regular solutions 
can be found and the system capsizes into chaos. Thus 
these "nonoptimal" solutions seem to be unstable. 

To analyze the stability of the found solutions we 
present the N$(\) dependence in Fig.[T](b), for d = 2, and 
in Fig. [5{b), for d = 3. Note these curves were built for 
"guess" functions corresponding to a minimal A. Apply- 
ing the Vakhitov-Kolokolov criterion 17] to the results 
shown on these plots one can conclude that the presented 
solutions are stable in 2D case, whereas some unstable 
solitons can exist in three dimensions. For d — 3, a sta- 
ble solution can be generated starting from a threshold 
plasmon number A cr s» 1100 and at |A| > 0.05 (see also 
the discussion in Sec. IIV[) . We show an example of a 
unstable soliton in Fig. [3J 

In Fig. Hlb), f° r d = 2, one can see that the new 
term in Eq. (BJ) , - (d - l)/x x tpd\ij;\ 2 / dx, does not 
produce any significant effect (compare solid and dash- 
dotted lines). In 3D case, Fig.[2jb), the difference is just 
quantitative: the critical plasmon number and critical 
frequency are shifted. Therefore our results are in agree- 
ment with Ref. fnl where NLSE with a nonlocal term 
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FIG. 1. The analysis of solutions of Eq. (12 . 23 fl in 2D case, (a) The function A(cr) obtained on the basis of Eq. (|2.23[) 
using a "guess" function ip = ip g (|3.1[) for different values of TVo. (b) The dependence N^,(X): solid line corresponds to a 
numerical solution of Eq. (|2,23[) . dashed line represents an "analytical" curve corresponding to minimal A a , and dash-dotted 
line corresponds to a numerical solution of Eq. (|2.23|l at the absence of the term (d — l)/x x ijjd\ip\ 2 /dx (see Ref. [HI]), 
(c) Examples of solutions of Eq. (|2.23[) obtained for No = 100; solid line: ~ 113, A ~ —0.20, and a ~ 3.83 (optimal 
parameters minimizing A), dashed line: m 70, A ~ —0.06, and a ~ 9.56, and dash-dotted line: « 595, A ~ —0.06, and 
a ~ 2.83. (d) The same as in panel (c) for No = 200; solid line: N^, « 210, A ~ —0.55, and a ~ 3.55 (optimal parameters), 
dashed line: « 113, A -0.20, and a w 7.91, and dash-dotted: w 113, A w -0.20, and a « 2.64. Note that in the 
latter case dashed and dash-dotted lines practically coincide. 



was analyzed. One should also notice that the numeri- 
cal curve in Fig. [TJb) (solid line) is in a good agreement 
with the "analytical" X(N^) dependence (dashed line), 
A a = -(AV - 87r) 2 /(887rAV), found from Eq. (l2~23|) with 
help of the "guess" function (|3.1[) . 



IV. APPLICATIONS 

Stable spatial solitons, involving nonlocal nonlineari- 
ties, similar to plasma structures described in the present 
work, were reported to be obtained in various labora- 



tory experiments in plasma physics, condescended mat- 
ter, and nonlinear optics (see, e.g., reviews 043 ^md 
references therein). Another opportunity for the phys- 
ical realization of plasmoids described in Sees. [TT1 and Mil 
consists in their implementation as a rare natural at- 
mospheric electricity phenomenon called a ball lightning 
(BL) 0. 

According to the BL observations, most frequently it 
has a form of a rather regular sphere with the diame- 
ter of (20 — 50) cm [l9j]. However big BLs with the size 
more than one meter were reported to exist (20| . Besides 
spherical BL, snake-like objects were observed 20]. The 
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(c) (d) 

FIG. 2. Panels (a) and (b) are the same as in Fig. [Ua,b) but correspond to 3D case, (c) Examples of stable solutions of 
Eq. (f!T2"3j) [right-handed branch in panel (b)] obtained for No = 1100; solid line: « 1150, A « -0.10, and a » 3.90 (optimal 
parameters), dashed line: N$ ~ 1429, A ~ —0.01, and a ~ 8.22, and dash-dotted line: N$ ~ 1219, A ~ —0.01, and a « 3.10. 
(d) The same as in panel (c) for No = 1200; solid line: N$ « 1197, A ~ —0.13, and cr ~ 3.80 (optimal parameters), dashed line: 
N^ » 1429, A » -0.01, and cr « 8.82, and dash-dotted line: iV,/, « 4354, A « -0.01, and a « 3.04. 



lifetime of BL can be up to several minutes [T3| . The es- 
timates of energy of BL were obtained only in cases when 
it produced some destruction while disappearing. These 
estimates give for the energy of BL the value in the range 
(several kJ - several MJ) [21(. However, since in many 
cases BL disappears just fading, one should assume that 
a very low energy BL can exist. 

Numerous BL models are reviewed in Ref. [22| . Despite 
very exotic theoretical descriptions of BL were proposed, 
it is most probable that this object is a plasma based 
phenomenon. In Refs. [23, l24[ we developed a BL model 
based on radial oscillations of electron gas in plasma, 
studying these oscillations in both classical and quantum 
approaches. The present work is a development of our BL 
model. As in our previous studies, here we also treat BL 
on the basis of radial oscillations of electrons. However, 



as we show in Sees. |TT] and IIII1 local and nonlocal elec- 
tron nonlinearities as well as the interaction of electrons 
with ions play an important role for the stability of an 
axially or spherically symmetric plasmoid. For example, 
the inhomogeneity of ion density was not accounted for 
in Ref. [23§ ■ It should be noted that the idea that spher- 
ically symmetric oscillations of electrons underlie the BL 
phenomenon was independently put forward in Ref. [25j j . 

In Fig. 2] we present the characteristics of axially sym- 
metric (snake- like BL [20]) and spherically symmetric 
(spherical BL) plasmoids calculated on the basis of the 
results of Sec. lIIII The total energy of electric field inside 
a plasmoid (energy density in 2D case) and the effective 
plasmoid radius are defined as 

W E = ljE^V, Rln = ^-Jr^V. (4.1) 
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X 

FIG. 3. An example of the unstable solution of Eq. (|2.23[) 
[left-handed branch in Fig. Hfb)] obtained for No = 770. It 
corresponds to N$ « 1174, A w -0.02, and a « 4.20. 

The frequency of the total oscillation, including the main 
harmonic oscillation which was separated in the deriva- 
tion of Eq. (gHHJ), is / = f p - |A/|, where f p = 2.8 x 
Y^MIlz is the Langmuir frequency. Here we take into 
account that frequency shift A/ for a stable plasmoid is 
negative [see Figs.QJa) and^a)]. 

One can see in Fig. 2] that in frames of our BL model 
we predict the existence of low energy (We > erg/cm in 
2D case and We > 10 2 erg in 3D case) plasma structures 
with the typical size R c ft < 1.5 m. This kind of plasmoids 
should appear in low density n e ~ 10 5 cm~ 3 and hot 
T e ~ 10 6 K plasma. Such a density of plasma can well 
exist in the Earth ionosphere [26j . A linear lightning 
can provide the plasma heating up to T e ~ 10 6 K since, 
e.g., Ti ~ 10 4 K in the lightning channel 27]. Thus we 
can assume that such high temperatures of plasma can 
be present in some localized area of ionosphere during a 
thunderstorm, making possible the existence of plasma 
structures described in the present work. It is interesting 
to notice that several meters objects were reported to 
appear in the BL observations near airplanes [28[. 

To demonstrate the stability of the described plas- 
moids with respect to decay, in Fig. [5] we present the ratio 
W E /W k versus W E , where W k = (3/2)n T e V cS is the to- 
tal thermal energy of a plasmoid (energy density in 2D 
case) and V e s is the effective volume of a plasmoid, which 
is equal to 7ri?g ff for d = 2 and to (4/3)7r_R^ ff for d = 3. 
Indeed, hot electrons with T e ~ 10 6 K could just escape 
the plasmoid volume making it unstable. One can how- 
ever see in Fig. [5] that for stable plasma structures the 
ratio Wft/Wk has a tendency to increase reaching the 
unit value at a certain soliton energy (We ~ 1 erg/cm 
for d = 2 and We ~ 215 erg for d = 3). It means 
that hot electrons will not escape form the plasmoid vol- 
ume as soon its energy has this critical value. It should 
be noticed that an unstable plasmoid, corresponding to 
the lower branch in Fig. [SJb), will lose electrons since 



We/Wc < 1 always. 

Comparing the predicted plasmoid properties with the 
characteristics of circumterrestrial (not ionospheric) BL 
one can say that our model is not directly applicable for 
the description of such an object since its energy > 10 kJ 
and the size ~ (20 — 50) cm are beyond our predictions. 
We can however use our model at the initial stages of the 
plasmoid formation when the amplitude of the electric 
field is so high since in Eq. (j2 . 18[) we keep only cubic 
nonlincarity. Higher nonlinear terms, which seem to be 
important for denser plasma, can explain smaller size and 
bigger energy content. 

A natural plasmoid appearing in circumterrestrial at- 
mosphere is surrounded by the neutral gas. It means that 
plasma in the interior of a plasmoid should be maintained 
in the state with proper ionization during its lifetime. It 
is however known [29| that plasma of a low energy plas- 
moid, without an internal energy source, will lose energy 
and recombine back to a neutral gas in the millisecond 
time scale at the atmospheric pressure. Thus the lifetime 
of such an object will be extremely short. It was sug- 
gested that under certain conditions plasma can reveal 
superconducting properties [301 ] . This mechanism could 
prevent the energy losses and thus the recombination, 
providing the long lifetime of a low energy plasmoid. 



V. CONCLUSION 

In conclusion we mention that in the present work we 
have studied axially and spherically symmetric Lang- 
miur solitons in warm plasma. In Sec. |TT] we started 
with the discussion of radially symmetric oscillations of 
electrons in plasma and then separated the motions on 
different time scales. In Sec. Ill Al we briefly reviewed 
the previous works on the theory of stable plasma struc- 
tures which involved local electron nonlinearities. Then, 
in Sec. Ill Bl we have calculated the contribution of the 
slowly varying electron density to the electron pressure 
and derived, together with the results of Ref. [ll(, the 
new NLSE (|2.18|) which accounts for nonlocal electron 
nonlineariries. These additional nonlinear terms do not 
disappear in one dimensional case. This fact was also 
demonstrated in Appendix lAl using Lagrange variables. 

The solutions of this new NLSE, rewritten in the equiv- 
alent form (|2.23[) for the radially symmetric case, have 
been analyzed numerically in Sec. IIIII We presented 
the examples of some of the solutions of Eq. (|2.23p [see 
Fig. [Hc,d) and Fig. [Hc,d)] and analyzed their stability. 
It has been found that solutions in 2D case are stable, 
whereas in tree dimensions some unstable solitons can ex- 
ist. We also compared out results with Ref. [ll|, where 
analogous NLSE was considered. 

In Sec.|TV]we suggested that the described axially and 
spherically symmetric solitons can be realized in the form 
of a natural plasmoid, a ball lightning. We have com- 
puted the characteristics of such a plasma structure, like 
energy and radius (see Fig. 2]), predicted in frames of our 
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FIG. 4. Physical characteristics of axial [panels (a) and (b)] and spherical [panels (c) and (d)] plasmoids expressed in terms 
of Te = T e /10 6 K and n 5 = no/10 5 cm -3 , (a) The energy density of electric field vs. the frequency shift, (b) The effective 
radius vs. the energy density, (c) The total energy of electric field vs. the frequency shift. Stable plasmoids correspond to the 
right-handed branch, (d) The effective radius vs. the total energy. Stable plasmoids correspond to the lower branch. 




FIG. 5. The ratio Ws/Wk vs. We for (a) 2D plasmoid and (b) 3D plasmoid. In the panel (b) a stable plasmoid corresponds 
to the upper branch and an unstable one to the lower branch. 
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model. Comparing these characteristics with the proper- 
ties of a natural plasmoid one concludes that our model 
describes a low energy BL which can exist in hot iono- 
spheric plasma. Note that in many cases BLs were ob- 
served form airplanes |28| at high altitudes. We also ex- 
amined the stability of plasmoids with respect to decay 
because of the escape of hot electrons. 

We have already considered a BL model based on 
radial oscillations of electrons in our previous publica- 
tions [23|,[24]]. In the present work we performed more de- 
tailed analysis of electron nonlinearities and showed that 
they play a crucial role for the stability of a plasmoid. Al- 
though plasma structures described in the present work 
does not reproduce all the properties of circumterrestrial 
(not ionospheric) BL [lj| , we can use our results for 
the description of the BL formation, when the amplitude 
of the electric field is not so big. 

As we have found in Sec. IIIII [see also Fig. Eta,c)], the 
total frequency of electron oscillations is less than Lang- 
muir frequency, since the frequency shift A/ is negative. 
It is important fact for the future experimental studies of 
BL. For example, the electron density in a linear lightning 
discharge can be ~ 10 17 cm -3 [3l|, giving for the plasma 
frequency a huge value of ~ 10 3 GHz. In a laboratory it 
is extremely difficult to create a strong electromagnetic 
field of such a frequency to generate BL. However, if the 
frequency of electron oscillations has a tendency to de- 
crease, it can facilitate the plasmoid generation. 
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Appendix A: Lagrange variables 

The propagation of waves in plasma was analyzed in 
Sec. Ill Al using Euler variables which are more convenient 
for the practical purposes. There is, however, a Lagrange 
approach for the treatment of plasma waves. Instead of 
describing of plasma characteristics, like velocity, den- 
sity etc, in a certain point of space, one can study the 
dependence of these quantities on the initial coordinates 
of plasma particles. 

Let us change the variables (r, t) —> (p, r) in Eqs. (12 . If) 
and (|2.5p as 



d_ 

dp 



dr d d 
dp dr ' dr 



d_ 
di 



d_ 

' dr 1 



(Al) 



where p is the initial coordinate of an electron, £ is de- 
viation of an electron from the equilibrium, and r is the 



new temporal variable. An electron is supposed to be in 
the equilibrium initially, £(p, 0) = 0. In Eq. (|A1|) we use 
the definition of electron velocity, dr/dr = d^/dr = v e . 
Unlike the Euler picture, the continuity equation in La- 
grange variables can be integrated and it does not contain 
the time derivative, 



noP 



n e r d 

dp 



(A2) 



where tiq = n e (r = 0) is the initial (unperturbed) elec- 
tron density, which is supposed to be uniform. 

Using Eqs. (|A"Tj) and (|A"2"j) we can obtain from Eq. (j2~T) 
a single nonlinear equation for electron velocity, 



n 



2r"v' 



d-l 

>e 

r 



(— 


-;)] 







0. 



(A3) 



where a "dot" and a "prime" mean the derivatives with 
respect to r and p. To derive Eq. (|A3|) we assume that 
electrons obey the adiabatic equation p/ril" — const, 
where p is electron pressure [pij = pSij, see Eq. (12.21) ] and 
7e = 3 since we study electrostatic oscillations (32( • Anal- 
ogous assumptions were taken in Ref. [l3j] to study the 
nonlinear one-dimensional plasma waves in warm plasma 
within the Lagrange picture. 

Note that Eq. (|A3j) is an exact one, which does not 
suppose any expansion over a small parameter. We can, 
however, discuss the small deviations from the equilib- 
rium position, 



i\ i 1 + — ) o e — 3v^-7r- 



n 



dp 



1 d 



pd- 1 dp 



d-l 



■ v e v e - 



p 



3i4 



+(d-i)( d -z) v 4 + (d-i) f 2 ^-^ 

p A V p p 



p 



0, 



(A4) 



keeping only quadratic nonlinearities. As in Sec. Ill Al 
n — rii — no, stays for the perturbation of the ion density 
in Eq. ([A"4]) . 

It can be noticed that Eq. (|A4|) has analogous structure 
as Eq. (|2.18l) before the separation of the main harmonic 
(for the details see, e.g., Ref. [Hj]). We can see, however, 
that nonlinear terms do not disappear completely at d = 
1 if we consider warm plasma with T e ^= 0. One can also 
notice that these nonvanishing nonlinear terms arise from 
the pressure term in the initial plasma hydrodynamics 
equations (I2.1[) . Analogous conclusion was obtained in 
Sec. lIIBI using Euler coordinates. Note that the nonlinear 
plasma waves were also studied in Refs. [Ill, [34]]. It was 
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found in Ref. [13jj that the nonlinear terms are important 
even for one dimensional plasma oscillations for the case 
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